##use package in huberโs m estimate library(robustreg) zzxxc<-robustRegH(datatata\(value ~ datatata\)x1+datatata$x2, data=datatata,tune=1.345,m=TRUE,max.it=1000,tol=1e-5 ) zzxxc[[1]]